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Nucleon structure study is one of the most important research areas in modern physics 
and has challenged us for decades. Spin has played an essential role and often brought 
surprises and puzzles to the investigation of the nucleon structure and the strong interac- 
tion. New experimental data on nucleon spin structure at low to intermediate momentum 
transfers combined with existing high momentum transfer data offer a comprehensive 
picture in the strong region of the interaction and of the transition region from the strong 
to the asymptotic-free region. Insight for some aspects of the theory for the strong in- 
teraction, Quantum Chromodynamics (QCD), is gained by exploring lower moments of 
spin structure functions and their corresponding sum rules (i.e. the Bjorken, Burkhardt- 
Cottingham, Gerasimov-Drell-Hearn (GDH), and the generalized GDH). These moments 
are expressed in terms of an operator-product expansion using quark and gluon degrees 
of freedom at moderately large momentum transfers. The higher-twist contributions 
have been examined through the evolution of these moments as the momentum transfer 
varies from higher to lower values. Furthermore, QCD-inspired low-energy effective the- 
ories, which explicitly include chiral symmetry breaking, are tested at low momentum 
transfers. The validity of these theories is further examined as the momentum transfer 
increases to moderate values. It is found that chiral perturbation theory calculations 
agree reasonably well with the first moment of the spin structure function gi at low 
momentum transfer of 0.05 - 0.1 GeV^ but fail to reproduce some of the higher mo- 
ments, noticeably, the neutron data in the case of the generalized polarizability Sj^x- 
The Burkhardt-Cottingham sum rule has been verified with good accuracy in a wide 
range of assuming that no singular behavior of the structure functions is present at 
very high excitation energies. 
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1. Introduction 

The nucleons (proton and neutron), the basic building blocks of nuclear matter, 
account for more than 99.9% of the mass of the visible matter in the universe. 
Understanding the nucleon structure is one of the most important issues in modern 
physics and has challenged us for decades. The internal structure of the nucleon 
presents a wealth of fundamental questions, that have a deep impact on our under- 
standing of Nature and the universe we live in. 

1 
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The strong interaction, that is mostly responsible for the structure of the nu- 
cleon, poses a special challenge to us due to the strong and non-linear nature of 
the force. Quantum Chromodynamics (QCD) is the accepted theory for the strong 
interaction. Quarks and gluons are the fundamental particles and force mediators, 
while nucleons and mesons (pions, kaons, etc.), collectively called hadrons, are com- 
posite particles and effective force mediators. QCD has been well tested at high 
energy (asymptotic-free) when the interaction is relatively weak due to the running 
of strong coupling constant. However, at low energy (scale of the size of a hadron), 
the strong force becomes truly strong, leading to chiral symmetry breaking and 
the confinement of the quarks and gluons inside the hadron where none can escape 
to become free particles. In the strong region of QCD, the usual calculation tech- 
nique of perturbative expansion fails, and QCD can not be easily solved. Effective 

rn 

field theories, such as chiral perturbation theory are often used. Recently, great 
progress has been made in the advancement of lattice QCD which has begun 
to provide some calculations of QCD in the strong region. Another approach^ is 
to stay with the continuous field theory by using the Schwinger-Dyson equations. 
Newly developed anti-de Sitter /conformal field theory (AdS/CFT) ''^ might provide 
another useful tool to study strong QCD. 

In addition to the truly strong region, understanding the transition from weak 
(perturbative) to strong (non-perturbative) region is also very important. One ap- 
proach is to start from high energy where the quarks and gluons (collectively called 
partons) ar e asym ptotically free. With the help of the Operator Product Expan- 
sion (OPE)ISlSlxl, one can systematically expand to a twist series, where the leading 
twist term gives the parton distributions and the higher-twist terms characterize 
the quark-gluon and quark-quark correlations. Lattice QCD can also provide use- 
ful calculations. Constituent-quark models are also often used to provide intuitive 
pictures in this region. 

Experimentally, deep-inelastic lepton scattering (DIS) together with other pro- 
cesses, such as Drell-Yan, direct photon and heavy quark productions, have been 
powerful tools to study nucleon structure. The unpolarized structure functions ex- 
tracted from the last five decades, covering five-order of magnitudes in both Q^ — 
and x— ranges, are the most extensive set of precision data in nucleon structure 
study. The longitudinally polarized structure functions extracted from the last three 
decades have reached reasonable precision, covering two to three-order of magni- 
tudes in Q^— and x— ranges. Transversely polarized measurements have been the 
recent focus, but are still scarce. The lower moments of the structure functions pro- 
vide the most direct tests and comparisons with theoretical calculations through 
sum rules. 

In the study of nucleon structure, spin has played an essential role and often 
brought big surprises. Spin was first introduced as an internal property of a particle, 
related to the magnetic moment, to explain the famous Stern-Gerlach results 
that the silver atoms split into two beams after passing through an inhomogeneous 
magnetic field. The discovery of the anomalous magnetic moment of the proton^ 
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was the first big spin surprise. This discovery was the beginning of the (indirect) 
study of nucleon structure, since the anomalous magnetic moment was one piece of 
evidence for the internal structure of the nucleon. Later on, it was related to the 
integration of the excitation spectrur n. of nu cleon spin structure by the sum rule of 
Gerasimov, Drell and Hearn (GDH) The direct study of nucleon structure 

started decades later when Hofstadter^I^used elastic electron scattering to directly 
measure the form factors of the proton. 

In the last thirty years, the spin structure of the nucleon has led to very 
productive experimental and theoretical activity with exciting results and new 
challenge This investigation has included a variety of aspects such as testing 
QCD in its perturbative region via spin sum rules (like the Bjorken sum rule^^ 
and understanding how the spin of the nucleon is built from the intrinsic degrees 
of freedom of the theory: quarks and gluons. Recently, results from a new gener- 
ation of experiments performed at Jefferson Lab seeking to probe the theory in 
its non-perturbative and transition regions have reached a mature state. The low 
momentum-transfer results offer insight in a region characterized by the collective 
behavior of the nucleon constituents and their interactions. In this region, it has 
been more economical to describe the nucleon using effective degrees of freedom like 
mesons and constituent quarks rather than current quarks and gluons. Furthermore, 
distinct features seen in the nucleon response to the electromagnetic probe, depend- 
ing on the resolution of the probe, point clearly to different regions of description, 
i.e. a scaling region where quark-gluon correlations are suppressed versus a coher- 
ent region where long-range interactions give rise to the static properties of the 
nucleon. 

In this review we describe an investigation ttSHSU of the spin structure of the 
nucleon through measurements of the helicity-dependent photoabsorption cross sec- 
tions or asymmetries by using virtual photons covering a wide resolution spectrum. 
These observables are used to extract the spin structure functions gi and 52 and 
evaluate their moments. These moments are key ingredients to test QCD sum rules 
and unravel some aspects of the quark-gluon structure of the nucleon. 



2. Formalism 

Leptons do not involve the strong interaction directly. With the electromagnetic in- 
teraction well-understood, lepton beams have been a very powerful precision probe 
to study nucleon structure. Consider deep-inelastic scattering of polarized leptons 
on polarized nucleons. We denote by m the lepton mass, k {k') the initial (final) 
lepton four-momentum and s {s') its covariant spin four-vector, such that s ■ k = 
(s' ■ k' = 0) and s ■ s = —1 {s' ■ s' = —1); the nucleon mass is M and the nucleon 
four-momentum and spin four-vector are, respectively, P and S. Assuming, as is 
usually done, one- photon exchange (see Fig. [T]), the differential cross section for 
detecting the final polarized lepton within the solid angle dfl and the final energy 
range (£", E' + dE') in the laboratory frame P = {M, 0), k = [E, k), k' = {E' , k'), 
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can be written as 



-— / 1/1/ A-"' 

dndE' 2Mq4 E 
where q = k — k' and a is the fine structure constant. 



lepton 



nucleon 




Fig. 1. Feynman diagram for deep-inelastic lepton-hadron scattering 



The leptonic tensor L^i, is given by 

L^^k, s; k') = [u{k', s') 7^ u{k, s)]* [u{k' , s') 7. u{k, s)] (2) 

s' 

and can be spht into symmetric {S) and antisymmetric {A) parts under fj, and v 
interchange: 

L^Ak,s-k')=2{L^J {k;k') + zL(^^ {k,s;k')} (3) 
where u{k,s) and u{k',s') arc electron spinors, and 

L^fJ {k; k') = kX + - ffM- • fc' - "i") (4) 
LI,-^) (/c, s; fc') = m (5) 

The unknown hadronic tensor W^i, describes the interaction between the virtual 
photon and the nucleon and depends upon four scalar structure functions: the 
unpolarizcd functions Fi^2 and the spin-dependent functions gi.2 (ignoring parity- 
violating interactions). These functions contain information on the internal struc- 
ture of the nucleon. They can be experimentally measured and then be studied in 
theoretical models, such as the QCD-parton model. The functions can only depend 
on the scalars and q ■ P. Usually people work with 

Q2 = -q^ and x = Q^/2q ■ P = Q'^/2Mu (6) 

where v = E — E' \s the energy of the virtual photon in the Lab frame. The variable 
X, known as "x-Bjorken" , is the fraction of momentum carried by the struck quark 
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in the simple parton model. We also refer to the invariant mass of the (unobserved) 
final state, W ^ {P + q)^ ^ + 2Mv - Q"^. 
Analogous to Eq. ([3)) one has 



W^^.(g; P, S) = Wll\q- P) + I Wlt\q; P, S). 
The symmetric part, relevant to unpolarized DIS, is given by 



Wj,t\q-P)^2 



2 



(7) 
(8) 



P.q 



p.- 



P.q 



The antisymmetric part, relevant for polarized DIS, can be written as 
Wlt\q- P, s) = 2 e^.^p q4 NPsfGi{v, Q^) + [MvS^ - (S ■ q) P^] G^i^, Q^) 



or, in terms of the scaling functions 171.2 



(9) 

(10) 



wl,t\q\P,'- 



2M 
'P~q 



■ . {S-q)P^ 



52(x,g2) . 



Note that these expressions are electromagnetically gauge-invariant: 



qf'W, 



0. 



In the Bjorken limit, or deep-inelastic regime. 



— )■ 00 



— !> 00 , X fixed. 



(12) 



(13) 



the structure functions Pi_2 and gi^2 are known to approximately scale, i.e., vary 
very slowly with at fixed x - in the simple parton model they scale exactly; in 
QCD their Q^-evolution can be calculated perturbatively. 

The differential cross section for unpolarized leptons scattering off an unpolar- 
ized target can be written as 



-{k,s,P,^S;k') 



2A/q4 E 



dn dE' 

While differences of cross sections with opposite target spins are given by 



(14) 



dVL dE 



-{Ks,P,~S-k') 



d^a 



dn dE 



- {k,s,P,S;k') 



(15) 



After some algebra one obtains from Eqs. (|4I5I8I11I14I15P expressions for 
the unpolarized cross section and polarized cross-section differences (note that the 
lepton-mass terms are neglected): 
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—F,{x,Q')ta.n'- + -F2{x,Q') (16) 



The unpolarized cross section is given by 

dO, dE' \dQ,J Mott 
where 9 is the scattering angle in the laboratory frame and the Mott cross 
section is given by 

/daN ^ a" cos^ f 

\dnJMott 4£2sin^|' ^ 

For the lepton and target nucleon polarized longitudinally, i.e. along or op- 
posite to the direction of the lepton beam, the cross-section difference under 
reversal of the nucleon's spin direction (indicated by the double arrow) is 
given by 



rf2<^ d''^ WE' 



dQ dE' dQ dE' uEQ^ 



(E + E'cose) gi{x,Q^)-2Mxg2{x,Q^) . 

(18) 



• For nucleons polarized transversely in the scattering plane, one finds 
dV^-f d^cT^JJ- Aa'^E''^ r , ^ME , ^2^1 

rfjyrf^-^orf^ = -W K'^ ^ + ^^^^"'^ )J • ^^^^ 

These two independent observables allow measurement of both g\ and 52 (as has 
been done at SLAC and in Jefferson Lab's Halls A and C), but the transverse cross 
section difference is generally smaller because of kinematic factors and therefore 
more difficult to measure. Only in the past decade has it been possible to gather 
precise information on g^-, which indeed turns out to be usually smaller than g\ in 
the deep-inelastic region. 

Experimental results are often presented in the form of asymmetries, which are 
ratios of the cross-section differences to the unpolarized cross section. 

For a longitudinally polarized target, the measured asymmetry is 

A ~ ~ (20) 

2 dc unpaid 



and for a transversely polarized target 

da^^ - da^^ 



(21) 



It is customary to introduce the (virtual) photon-nucleon asymmetries A\^2' 

= ^ = ^1^^ (22) 

and 

A2 = ^=7[^], (23) 
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where 



7 



2Mx Q 



(24) 



Q V ' 

and tJT, ott and oj^t are the virtual photon cross sections (see below for the 
definitions of virtual photon cross sections). 
From Eqs. (116118119^ . it foUows that 



and 



where 



D 



All =D{Ax+T]A2) 

y[(l + 7V2)(2-y)-2yW/Q^] 
j/2(l - 2m2/Q2)(l + 72) + 2(1 + i?)(l -y- -f'^y'^/A) 

[1 + 7'y/2(l + 2m2y/Q2)]^i_ 2^.^2^2/4- 



7? = 7 



(l-y/2)(l + 72y/2)-y2m2/Q2 

[i-y-y'(7V4 + ™VQ')] 



D 



(l-y/2)(l + 72t//2)-y2,„2/Q2 
1 - J//2 - y^m^lQ'^ 



l + -l'^y/2{l + 2m'^y/Q'^Y 



where 



and 



is the ratio of the longitudinal and transv erse vi rtual photon cross sections. 
There exists a restrictive bound on A2 l ^^ l ^^ l; 



IA2I < ^R{l + Ai)l2. 



(25) 
(26) 

(27) 

(28) 

(29) 
(30) 

(31) 
(32) 

(33) 



In the virtual photon notation, the inclusive inelastic cross section can be written 



in terms of a virtual photon flux factor and four partial cross section 



(Pa 

dn dE' 



(34) 



a = ax + ecTL - hPx\/2e{l - e)aLT - hPz\/T~~e^aTT, 



(35) 
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with the photon polarization 

^ " l + 2(l + i/2/Q2)tan2(0/2) ' ^^^^ 

and the flux factor 

aE'K 

- 2^W(l-e) ' ^^^^ 

where K is the " equ ivalent photon energy" and here we will use the definition 
according to Hand^ 

K = v{l-x) = —— . (38) 

h = ±1 refers to the two helicity states of the (relativistic) lepton. Pz is the target 
polarization along the direction of the virtual photon and Px , perpendicular to that 
direction in the scattering plane of the electron. 

The partial cross sections are related to the structure functions as follows: 

CTt = cri/2 + (73/2 = ' (39) 

(JTT = CTi/2 - CT3/2 = -j^yai -1 92) , (41) 

^LT = -J^^i9i + 92) , (42) 

where ai/2 and (73/2 are the helicity cross sections with 1/2 and 3/2 referring to 
the total helicity of the (virtual) photon-nucleon system. 

Note that in the above we have kept terms of order /Q^ and smaller. They 
are sometimes necessary in order to extract the correct experimental values of the 
structure functions from the measured asymmetries. However, the QCD analysis of 
the structure functions is often carried out at leading twist only, ignoring higher- 
twist terms, i.e., higher order in M^/Q^. This is clearly inconsistent in cases where 
the above terms could be important, for example Jefferson Lab experiments. 



3. Sum rules and Moments 

Sum rules involving the spin structure of the nucleon offer an important opportunity 
to study QCD. In recent years the Bjorken sum rule at large and the GDH sum 
rule at Q2 = have attracted large experimental'^'^ '^^ '^^ and theoreticaPS efforts 
that have provided us with rich information. This first type of sum rules relates 
the moments of the spin structure functions (or, equivalently, the spin-dependent 
virtual photon cross sections) to the nucleon's static properties. The second type 
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of s um rul es, such as the generahzed GDH sum ruld^lESl or the polarizabihty sum 
ruleJM, 

relate the moments of the spin structure functions to real or virtual 
Compton amplitudes, which can be calculated theoretically. Both types of sum 



rules are based on "unsubtracted" dispersion relations and the optical theorerrPSl^ 
The first type of sum rules uses one more general assumption such as a low-ene rgy 
theorenPlforthe GDH sum rule and the Operator Production Expansion (OPE/^ 
for the Bjorken sum rule to relate the Compt on amp litude to a static property. 

The formulation below follows closely Refs.SSEZl, Consider the forward doubly- 
virtual Compton scattering (VVCS) of a virtual photon with space-like four- 
momentum q, i.e., — ~ (f' — < 0. The absorption of a virtual photon on 
a nucleon is related to inclusive electron scattering. As discussed in the previous 
section, the inclusive cross section contains four partial cross sections (or structure 
functions): ctt, <Jl, cttt, o'lt, (or Fi, F2, gi, .92)- In this review, we will concentrate 
on the spin-dependent functions, cttt, ctlt (or gi, (72)- In the following discussion, 
we will start with the general situation, i.e., sum rules valid for all , then discuss 
the two limiting cases at low and at high Q^. 

Considering the spin-flip VVCS amplitude grr and assuming it has an appropri- 
ate convergence behavior at high energy, an unsubtracted dispersion relation leads 
to the following equation for grr'- 

p r I n^. pole, ^2^1 , ^ ^n. T :^£KQ>TlV^ , , 

where g^*^"^ is the nucleon pole (elastic) contribution, V denotes the principal value 
integral. The lower limit of the integration vq is the pion-production threshold on 
the nucleon. A low-energy expansion gives: 

Re[gTT{>y, Q') ~ g^rfi'^^ Q')] = {^)lTT{Q'y + ItAQ^V + 0{i^% (44) 

Itt{Q^) is the coefficient of the 0{i') term of the Compton amplitude. Eq. (|44p 
defines the g enerali zed forward spin polarizabihty 7rr(Q^) (or 7o(Q^) as it was 
used in Refs. 

USEi). Combining Eqs. ^ and the 0{ v) term yields a sum 

rule for the generalized GDH integral'^'^ 4'^: 

Itt{Q ) = -7-^ I dv 



2M^ 



1 fXO 







2 



gi{x, Q2) - ^M:^x-'g2{x, Q'^) dx. (45) 



As — ^ 0, the low-energy theorem relates 1(0) to the anomalous m agnetic moment 
of the nucleon, k, and Eq. (P5|) becomes the original GDH sum rulJl^mi 

m = r . (46) 
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The 0{v^) term yields a sum rule for the generahzed forward spm 
polarizabilitjBSmi: 



1tt{Q) = (^) 



dv 



Xo 



AM' 



dx. 



(47) 



Considering the longitudinal-transverse interference amplitude qlt, with the 
same assumptions, one obtains: 

M9lt{v. Q^) - glf{v, Q2)] ^ (|^)q/^^(q2) ^ QS^^^Q^y + 0{v^) (48) 

where the 0(1) term leads to a sum rule for Ilt{Q^), which is related to the ct^t 
integral over the excitation spectrum: 



2A'p r° 



dv 



giix,Q^)+g2{x,Q^) 



dx. 



(49) 



The 0{i'^) term leads to the generalized longitudinal-transverse polarizabilitjH^S^ 



SMQ^) = (i) 



dv 



Xo 



giix,Q^)+g2ix,Q^) 



dx. 



(50) 



Alternatively, we can consider the covariant spin-dependent VVCS amplitudes 
Si and ^2, which are related to the spin- flip amplitudes gxr and gLT- 

vM 



Si{v,Q' 



v^ + Q^ 



-9lt{v,Q 



v^ + Q^ 

The dispersion relation with the same assumptions leads to 



M 



where the 0(1) term leads to a sum rule for Ii{Q^)' 



2M2 



gi{x,Q'^)dx. 



(51) 



(52) 



(53) 



The 0{v^) term leads to the generalized gi polarizability: 



5lt 



2a 



x'^gi{x, Q^)dx 
Itt{Q') ~ hiQ') 



(54) 
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For 5*2, assuming a Regge behavior 'At v oo given by 5*2 — i^"^ with a2 < 
— 1, the unsubtracted dispersion relations for ^2 and vS2, without the elastic pole 
subtraction, lead to a "super-convergence relation" that is valid for any value of 

fi 

(72(a;,Q2)dx = 0, (55) 

which is the Burkhardt-Cottinghani (BC) sum ruleP. This expression can also be 
written as 



0Af2 j-xo 1 , 



(56) 



where Fp and Fd are the Pauli and Dirac form factors for elastic e-N scattering. 
The low-energy expansion of the dispersion relation leads to 

Re [(i.^2(^^,Q'))-Kf°'^(i',Q'))] 
= 2a/2(Q2) - ^(/tt(Q') - /i(Q'))^' + ^laAQ'y + 0{i^% (57) 

where the 0{v^) term gives the generalized g2 polarizability: 

, 2^ 167raM2 „ ^ „ 

792 (Q ) = — -QQ j X g2(x,Q )dx 

= SMQ') - iMQ') + {itt{Q') - h (Q^)) . (58) 

At high the OPE EEEESl for the VVCS amplitude leads to the twist ex- 
pansion: 



Jo A 



'^^(^^^ (59) 



.=2,4,... iQ'Y^-'^^' 

with the coefficients fir related to nucleon matrix elements of operators of twist 
< r. Here twist is defined as the mass dimension minus the spin of an operator, 
and the Hr are perturbative series in as, the running strong coupling constant. Note 
that the application of the OPE requires a summation over all hadronic final states 
including the ground state at a; = 1. 

The leading- twist (twist-2) component, /i2, is determined by matrix elements 
of the axial vector operator ip^p,j5ip summed over quark flavors, where are the 
quark field operators. It can be decomposed into flavor triplet (5^1), octet (as) and 
singlet (AS) axial charges, 

M2(Q') = {±^fJA + ^as) + ^AS + 0{as{Q^j), (60) 

where -l-(-) corresponds to proton (neutron), and the 0{as) terms involve the 
Q^— evolution due to the QCD radiative effects that can be calculated from pertur- 
bative QCD. The triplet axial charge is obtained from neutron /3-decay, while the 
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octet axial charge can be extracted from hyperon weak-decay matrix elements as- 
suming SU(3) flavor symmetry. Within the quark-parton model AS is the amount 
of nucleon spin carried by the quarks. Deep Inelastic Scattering (DIS) experiments 
at large have extracted this quantity through a global analysis of the world 



datatiS. 

Eqs. (|59| and (|60)) . at leading twist and with the assumptions of SU(3j flavor 
symmetry and an unpolarized strange sea, lead to the Ellis- Jaffe sum rule The 
difference between the proton and the neutron gives the flavor non-singlet term: 

r?(Q') - r?(Q2) = I + o{a,) + 0(1/Q2), (61) 
6 

which becomes the Bjorken sum rule at the — oo limit. 

If the nucleon mass were zero, Hr would contain only a twist-r operator. The 
non-zero nucleon mass induces contributions to fXr from lower-twist operators. The 
twist-4 term contains a twist-2 contri bution, 02 , and a twist-3 contribution, ^2, in 
addition to /2, the twist-4 componenlPESlMl 

M4 = (02 + 4d2 + 4/2) /9. (62) 

The twist-2 matrix element 02 is 

02 ^{A-p-pA} = 1 ^ e2 (P, 51^;, -f^f^iOnD^^^]?, S) , (63) 

where Cg is the electric charge of a quark with flavor q, D'^ are the covariant deriva- 
tives and the parentheses {• • • } denote symmetrization of indices. The matrix ele- 
ment a2 is related to the second moment of the twist-2 part of gi: 



a2(g')=2 / dxx' gi{x,Q^) . 
Jo 



(64) 



Taking Eq. ([64]) as the definition of 02, it is now generalized to any including 
twist-2 and higher- twist contributions. At low Q^, the inelastic part of a2 is related 
to , which is the generalized gi polarizability: 

Note that at large Q^, the elastic contribution is negligible and 02 becomes 02- 
The twist-3 component, d2, is defined by the matrix element*^ 

d.s^^pMp^} = i ^ {P, S\gi^, F^'^'' ^^H,\P, S) , (66) 
q 

where g is the QCD coupling constant, F^^" — {l/2)e^^'^°'^ Fap is the dual gluon-field 
strength tensor, Fap are the gluon field operators and the parentheses [• • • ] denote 
antisymmetrization of indices. This matrix element depends explicitly on gauge 
(gluon) fields. The gauge fields can be replaced by quark fields using the QCD 
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equation of motion Then, the matrix element can be related to the second 
moments of the twist-3 part of gi and 52 '■ 

rf2(Q') dx x''{2gi{x,Q^)+Zg2{x,Q^)^ 

= ^I'dx x' (g2(x, - gY'^ix, Q')) , (67) 

where g^^ is the twist-2 part of 32 as derived by Wandzura and Wilczell^ 

g^^ix, Q') = -g,{x, Q^) + f dy^-^^^ . (68) 

Jx V 

The definition of ^2 with Eq. (|67p is generalized to all Q^. At low Q^, the inelastic 
part of d2{Q^) is related to the polarizabilities: 

- ^^(7..+^7.J = ^{h{Q')-iTT{Q') + ^SMQ'))- (69) 

At large Q^, d2 becomes d2 since the elastic contribution becomes negligible. 
The twist-4 contribution to /i4 is defined by the matrix element 

/2 M^^^ = IT.4 ^\9^i ^) • (70) 

This matrix element depends also explicitly on gauge (gluon) fields. The gluon 
fields can be replac ed w ith the "bad" components of quark fields using the QCD 
equation of motion EESI, Then this matrix element can be related to moments of 
parton distributions: 

MQ^) = ^l dx x^{7gi{x,Q^) + 12g2{x,Q^)-9g3{x,Q^)y (71) 

where 

Jo g 

is the twist-4 distribution function defined in terms of the "bad" components of 
quark fields in the light-cone gauge, and p is along the direction of "bad" com- 
ponents (perpendicular to the light cone direction) 0. With only gi and g2 data 
available, /2 can be extracted through Eqs. (|59p and if the twist-6 or higher 
terms are not significant or are included in the extraction. 

The twist-3 and 4 operators describe the response of the collective color electric 
and magnetic fields to the spin of the nucleon. Expressing these matrix elements in 

^This 33 is not to be confused with the parity-violating structure function, which also uses the 
notion of ga in the literature. 
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terms of the components of F^'^ in the nucleon rest frame, one can r elate ^ 2 and 
/2 to color electric and magnetic polarizabihties. These are defined as^ ^^ l ^^ l 

XE 2M^S = {N\ ja X Ea \N) , XB 2M^S = (iV| jO Ba \N) , (72) 

where S is the nucleon spin vector, is the quark current, Ea and Ba are the 
color electric and magnetic fields, respectively. In terms of c?2 and /2 the color 
polarizabihties can be expressed as 

XE^l (2d2 + /2) , XB^l (4rf2 - /2) . (73) 

Recently, M. Burkardt pointed out that d2 and /2 are local correlators. 
He suggested that instead of calling them color polarizabihties, they may more 
appropriately be identified with the transverse component of the color-Lorentz force 
acting on the struck quark. Furthermore, he pointed out an interesting link between 
^2 and the Sivers function, a transverse-momentum-dependent (TMD) distribution 
function. Since this is beyond the scope of this article, interested readers can find 
the details in Ref. EH. 



4. Summary of experimental situation before JLab 

Before Jefferson Lab started experiments with polarized beams and targets, most of 
the spin structure mea sur ements of the nucleon we re performed at hig h- energy facil- 
ities like CERN (EMC EH and SMCl ^^ l ^9 | 60 | 61 | 62 |)^ DES Y (HERMES IS3IMI651) ^^^^ 
SLAC (ESOESI, EiaoEZl, EI42EHI, EMSISSI, E154l2niZl], E155 1221221 and E155x[22). 
The measured gi and 52 data were suitable for an analysis in terms of perturbative 
QCD. The impetus for performing these experiments on both the proton and the 
neutron was to test the Bjorkcn sum rule, a fundamental sum rule of QCD. After 
twenty-five years of active investigation this goal was accomplished with a test of 
this sum rule to better than 10%. The spin structure of the nucleon was unraveled 
in the same process. Among the highlights of this effort is the determination of 
the total spin content of the nucleon due to quarks, AE (see Eq. (|60|) ). This study 
also revealed the important role of the quark orbital angular momentum and gluon 
total angular momentum. The main results from this inclusive double spin asym- 
metry measurement program have led to new directions: namely the quest for an 
experimental determination of the orbital angular momentum contribution I 7'7 | 78 | 
(e.g., with Deep Virtual Compton Scattering and Transverse Target Single Spin 
Asymmetry measurements at Jefferson Lab, HERA and CERN) and the gluon spin 
contribution (with COMPASS and RHlC-spin I^Sl experiments). These efforts 
will be ongoing for the next few decades. 
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5. Recent results from Jefferson Lab 
5.1. J Lab experiments 

With a high current, high polarization electron beam of energy up to 6 GeV and 
state-of-the-art polarized targets, JLab has completed a number of experiments, 
which have extended the database on spin structure functions significantly both in 
kinem atic range (low and high x) and in precision. The neutron results are from 
Hall using polarized '^hJ^ as an effective polarized neutron target and two 
high resolution spectrometers. The polarized luminosity reached 10'^^ s~^cm~^ with 
in-beam polarization improved from 35% (1998) to over 65% (2008). The proton 
and deuteron results are from Hall eISII with the CLAS detector and Hall C with 
the HMS spectrometer and using polarized NH3 and ND3 targetJ^ with in-beam 
polarization of about 80% and 40%, respectively. 



5.2. Results of the generalized GDH sum for the neutron and ^He 

The spin structure functions gi and 172 (or axT and u^t) were measured in Hall A 
experiment E94-010 on ■^He from break-up threshold to = 2 GeV covering the 
Q2_range of 0.1-0.9 GeV^. The generalized GDH integrals I{Q^) (Eqgg) (open 
symbols) were extracted for ^He (top plot of Fig. [5]) and for the neutron 
(bottom). The solid squares include an estimate of the unmeasured high-energy 
part. 

The •^He results rise with decreasing . Since the GDH sum rule at Q-^ = 
predicts a large negative value, a drastic turn around should happen at lower 
than 0.1 GeV^. A simple model using MAID plus quasielastic contributions esti- 
mated from a P WIA model indeed shows the expected turn around. The data 
at low should be a good testing ground for few-body chiral perturbation theory 
calculations when they are available. 

The neutron results indicate a significant yet smooth variation of I{Q^) to in- 
creasingly negative values as varies from 0.9 GeV'^ towards zero. The data are 
more negative than the MAID model calculation Since the calculation only in- 
cludes contributions to I{Q^) iov W <2 GeV, the model should be compared with 
the open squares. The GDH sum rule prediction, J(0) = — 232.8 /xb, is indicated 
along with extensions to > using two next-to-leading order chiral perturba- 
tion theory (^PT) calculations: one of them using the heavy baryon approximation 
(HBxPT) 95 (dotted line) and the other relativistic baryon xPT (RBxPT) (dot- 
dashed line) . Shown with a brown band is a RB^PT prediction including resonance 
effects^. The large unc ertaint y is due to the resonance parameters used. 

Improved calculationJ^HEHI g^g .^^ell as further measurements will help clarify 
the situation. 
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Fig. 2. Results for the GDH sum /(Q^) for ^He (top) and the neutron (bottom). The ^He GDH 
results are compared with the MAID model plu s a guasielastic contributi on. The neutron GDH 
results are compared with xPT calcul atio ns of ref.lSSl (dotted curve) and ref . I2ZI (dot-dashed curve}. 
The MAID model calculation of ref. I^^J jg represented by a solid curve. Data from HERMES 
are also shown. 
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5.3. First moments of gi and the Bjorken sum 

The first moment of gi, Fi, was extracted from E94-010 for the neutron ^'^ and 
^He 11^, and from the CLAS egl experiment for the proton ^"^j the deuteron 1^ 
and the neutron (from the deuteron with the proton contribution subtracted) over 
a Q^— range from 0.05 to 5 GeV^. The results are plotted in Fig. |31 Also plotted 
are the results at = 1.3 GeV^ from the Hall C RSS experiment These 
moments show a strong yet smooth variation in the transition region (Q^ from 1 
to 0.05 GeV^). The lowest (0.05-0.1 GeV^) data are compared with two xPT 
calculations and are in reasonable agreement. 

Also plotted are the preliminary neutron data at very low of 0.04 to 0.24 
GeV^ from the Hall A E97-110 experiment. The precision of the preliminary data 
is currently limited by systematic uncertainties, which are expected to be reduced 
significantly when the final results become available. The new data are in good 
agreement with published results in the overlap region. The data in the very low 
Q^— region provide a benchmark test of the xPT calculations, since they are ex- 
pected to work in this region. The results agree well with both xPT calculations 
and also indicate a smooth transition as approaches ze ro. 

The bottom-right panel shows the first moment for p-nl221^ a flavor non-singlet 
combination, which is the Bjorken sum at large ^T\. The Bjorken sum was 
used to extract the strong coupling constant as at high (5 GeV^). An attempt 
was made to extract an "effective coupling" a^.g^ in the low Q^— region using the 
Bjorken sum, see section \5M 

At large Q^, the data agree with pQCD models where the higher- twist effects 
are small. 



5.4. First moment of g-^: B-C sum 

Measurements of g2 require transversely polarized targets. SLAG E155x I^H per- 
formed the only dedicated g2 measurement prior to JLab. At JLab, 52 and its 
moments have been extensively measured for the neutron with a polarized ' ^He 
target in a w ide range of k inematics in several Hall A experiments (E94-01ol^, 
E97-103[I^, E99-117ES1^ E02-012 1211 and E97-110 tlS). The measurement on 
the proton was performed in the RSS experiment in Hall C at an average of 
1.3 GeV2. 

The first moment of 52 is expected to be zero at all from the Burkhardt- 
Cottingham (B-C) sum rule. The SLAG E155x results yielded a first test of the 
B-G sum rule for the proton, deuteron and neutron (extracted from p and d data). 
The proton result (top panel of Fig. [4]) appears to be inconsistent with the B-G 
sum rule at the 2.75 a level. In addition to the large experimental uncertainty, 
there is an uncertainty associated with the low- a; extrapolation that is difficult to 
quantify. The deuteron and neutron (bottom panel of Fig. |4]) results from SLAG 
show agreement with the B-G sum rule with large uncertai nties. The most extensive 
measurements of the B-G sum rule come from experiments [ISHZHD qji ^jjg neutron 
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Fig. 3. Results of T^{Q'^) for p, d, n and p-n from the JLab Hall A, CLAS egl and Hall C RSS 
experiments. The slopes at = predicted by the GDH sum rule are given by the dotted line s. 
The dashed (plain) lines are the predictions from the Soffer-Teryaev IIQQI (Burkert-Ioffe 
model. The leading twist Q'^-evolution of the m ome nts is given by the gra v band. The inset shows 
comparisons with xPT calculations by Ji et al. and Bernard et al. 12ZI. 



using a longitudinally and transversely polarized '^He target in Hall A. The results 
for are plotted in the bottom panel of Fig. |4] in the measured region (open 
circles). The solid diamonds include the elastic contribution and an estimated DIS 
contribution assuming g2 = ■ The published results from E94-010 and the 

preliminary results from E97-110 and E02-012 are the most precise data on the B- 
C sum rule and are consistent with the expecta tion of zero, within small systematic 
uncertainties. The RSS experiment in Hall C took data on g2 for the proton 
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Fig. 4. The verification of the BC sum rule from Hall C RSS and Hall A experiments E94-010, E97- 
110 and E02-012 (neutron only), together with SLAC E155x data. The top plot is for the proton 
and the bottom for the neutron. The open symbols are the measured values and the solid symbols 
are the total moments, including the elastic and estimated contributions from the unmeasured 
high-energy region. 



and deuteron, at an average of about 1.3 GeV^. For the proton (top panel 
of Fig. [4|), the integral over the measured resonance region is negative, but after 
the elastic contribution and the estimated small- a; part are added, the preliminary 
result for the B-C sum is consistent with zero. Preliminary result on the deuteron 
and the neutron (from d and p) are also consistent with zero within uncertainties. 
The precision data from JLab are consistent with the B-C sum rule in all cases, 
indicating that §2 is a well-behaved function with good convergence at high energy 
though the high energy part is mostly unmeasured. 

5.5. Spin Polarizabilities: 70, 6lt 

The generalized spin polarizabilities provide benchmark tests of xPT calculations at 
low . Since the generalized polarizabilities have an extra l/v"^ weighting compared 
to the first moments (GDH sum or /lt), these integrals have less of a contribution 
from the large-i^ region and converge much faster, which minimizes the uncertainty 
due to the unmeasured region at large v. 

At lo w Q^, t he generalized polarizabilities have been evaluated with xPT cal- 
culations ISSHQSI. One issue in the xPT calculations is how to properly include the 
nucleon resonance contributions, especially the A resonance which usually domi- 
nates. As was pointed out in Ref. ^ while 70 is sensitive to resonances, 5lt is 
insensitive to the A resonance. Measurements of the generalized spin polarizabil- 
ities are an important step in understanding the dynamics of QCD in the chiral 
perturbation region. 
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The first results for tfie neutron generalized forw ard spin polarizabilities 7o(Q^) 
and Slt{Q^) were obtained at Jefferson Lab Hall A^ISlover a Q^— range from 0.1 
to 0.9 GeV^. We will focus on the low Q^— region where the comparison with xPT 
calculations is meaningful. 




Fig. 5. Left panels display results for the neutron spin polarizabilities 70 (top) and Sj^t (bottom). 
Solid squares are the results with statistical uncertainties. The dark bands indicates the systematic 
uncertainties. The dashed curves are t he H ByPT calculation^'-'''. The dot-dashed curves and the 
light bands are the RBxPT calculatioJSS! without and with'^''' the A and vector-meson contribu- 
tions, respectively. Solid curves are from the MAID model^^l The other panels are results for the 
spin polarizability 70 for the proton, p — n and p + n. 



The results for 7o(Q^) for the neutron are shown in the top-left panel of Fig. [5] 
for the two lowest Q^— values of 0.10 and 0.26 GeV^. The statistical uncertainties 
are smaller than the size of the symbols. The data are compared with a next- 
to-leading order, 0{p'^), HBxPT calculatioiJlQ^, a next-to-leading-order RB^PT 
calculatiorPSl and the same calculation explicitly including both the A resonance 
and vector-meson contributionJ23. Predictions from the MAID modeP^ are also 
shown. At the lowest point, the RBxPT calculation including the resonance con- 
tributions agrees with the experimental result. For the HB^PT calculation without 
explicit resonance contributions, the discrepancy is large even at — 0.1 GeV^. 
This might indicate the significance of the resonance contributions or a problem 
with the heavy baryon approximation at this Q^. The higher-Q^ data point is in 
good agreement with the MAID prediction, but the lowest data point at — 0.1 
GeV^ is significantly lower, consistent with what was observed for the generalized 
GDH integral result (Section 15. 2|) . Since Slt is insensitive to the dominating A- 
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resonance contribution, it was believed that 5lt should bemore suitable than 70 
to serve as a testing ground for the chiral dynamics of 

QQ]j 96 | 105 | rpj^g bottom-left 
panel of Fig. [5] shows 5lt compared to xPT calculations and the MAID predictions. 
It is surprising to see that the data are in significant disagreement with the xPT 
calculations even at the lowest-Q^ of 0.1 GeV^. This discrepancy presents a signifi- 
cant challenge to the present theoretical understanding. The MAID predictions are 
in good agreement with the results. 

The other panels of Fig. [5] present results of 70 for the proton ^22] ^nd for the 
isospin decompositions: p-n and p-|-nl31l. 

These results are in strong disagreements with both calculations. Since 

7o is sensitive to the A contributions, this may indicate that the treatment of the 
A contributions in the xPT calculations needs to be t aken into account properly. 
Progress in this direction is expected in the near future I^SI. 



5.6. Higher moment d-^ and higher twist: twist-3 

Another combination of the second moments, d2{Q^), provides an efficient way to 
study the high-Q^ behavior of nucleon spin structure, since at high it is related 
to a matrix element ^2 (color polarizability) which can be calculated from Lattice 
QCD. This moment also provides a means to study the transition from high to low 
Q^. In the left panel of Fig. HI d2{Q^) on the neutron is shown. The preliminary 
results from E02-012 ^211 and the published E94-010 results are shown as the 
solid squares and open circles, respectively. The bands represent the systematic 
uncertainties. The neutron resuhs from SLACp^open diamond) at high and 
from combined SLAG and JLab E99-117Eni(solid diamond) are also shown. The 
solid line is the MAID calculation containing only the resonance contribution. At 
moderate Q^, our data show that c?2 positive and decreases with Q^. At large 
q2 (q2 ^ 5 GeV^), data (SLAG alone and SLAG combined with JLab E99-117) 
are positive. The lattice QGD predictiorff^ at = 5 GcV^ is negative but close 
to zero and is about two sigmas away from the data. We note that all models 
(not shown at this scale) predict a small negative or zero value at large Q^. High- 
precision data at large will be crucial for a benchmark test of the lattice QGD 
predictions and for understanding the dynamics of the quar k-gl uon correlations. 

Also shown are the 6,2 result on the proton from Hall G at = 1.3 GeV^ 
and the same result evolved to = 5 GeV^ along with the SLAG data and the 
Lattice QGD calculation. 



5.7. Higher twist extractions: f-^, twist-4 

The higher-twist contributions to Fi can be obtained by a fit with an OPE series, 
Eq. (|59p . truncated to an order appropriate for the precision and Q^-span of the 
data. The goal is to determine the twist-4 matrix element f2- Once /Z4 is obtained, 
/2 is extracted by subtracting the leading-twist contributions of a2 and c?2 following 
Eq. ((62|l . The data precision and Q^— range allow a fit to both the /14 and /ig terms. 
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Fig. 6. The left plot shows the Hall A results of d2 for the neutron along with the world data at 
high Q'^, Lattice QCD, MAID model and HBxPT calculations. The right plot shows the Hall C 
result of d2 for the proton at = 1.3 GeV^, the same data evolved to = 5 GeV^ along with 
the SLAG data and lattice QCD calculation. 



To have an idea how the higher-twist terms (twist-8 and above) affect the twist-4 
term extraction, it is necessary to study the convergence of the expansion and to 
choose the Q^— range in a region where the fig, term is not significant. A fit with 
three terms (/i4, /ig and /is) was also preformed for this purpose. This study was 
possible only because of the availability of the high-precision low-Q^ data from 
JLab. 

Higher-twist analyses have been performed on the protorPH, the neutron^^ and 
the Bjorken surrP^l. An earlier proton analysis is availablJ^Hl but will not be pre- 
sented here, since that analysis used a different procedure. Fi at moderate is 
obtained as described in section 15.11 For consistency, the unmeasured low- a; parts 
of the JLab F^ and of the world data on Fi were re-evaluated using the same 
prescription previously used for F" and F^~". The elastic contribution, negligible 
above of 2 GeV^ but significant (especially fo r the proton) at lower values of 
Q^, was added using the parametrization of Ref. HQS. The leading- twist term fi2 
was determined by fitting the data at > 5 GeV^ assuming that higher twists are 
negligible in this Q^— region. A value of gA = 1.270 ± 0.045 was obtained for the 
Bjorken sum. Using the proton (neutron) data alone, and with input of qa from 
neutron beta decay and ag from hyperon decay (assuming SU(3) flavor symmetry), 
we obtained AE = 0.15 ±0.07 for the proton and AS = 0.35 ±0.08 for the neutron. 
We note that there is a significant difference between AS determined from the pro- 
ton and from the neutron data. This is the main reason why the extracted /14 and 
/2 from the Bjorken sum is different compared to the difference of those extracted 
individually from the proton and neutron, since the Bjorken sum does not need the 
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assumption of SU(3) flavor symmetry and AE was canceled. 

The fit results using an expansion up to (l/Q^) in determining /i4 are sum- 
marized in Table 1. In order to extract /2, as shown in Table 2, the target- 
mass correctio ns a -i were evaluated using the Blumlein-Boettcher world data 
parametrizatioiff^ for the proton and a fit to the world neu tron data, which in- 
cludes the recent high-precision neutron results at large aP^. The d2 values used 
are from SLAG E155>l^ (proton) and JLab E99-llfliS! (neutron). 



Table 1. Results of ^4, and /xg at = 1 GoV^ for proton, neutron and p~n. The uncertainties 
are first statistical, then systematic. 



Target 


Q'^ (GcV^) 




1,6 /M^ 




proton 
neutron 
p — n 


0.6-11.0 
0.5-11.0 
0.5-11.0 


-0.065±0.012 ± 0.048 
0.019±0.002± 0.024 
-0.060±0.045 ± 0.018 


0.143±0. 021 ±0.056 
-0.019±0.002± 0.017 
0.086±0.077 ± 0.032 


-0.026±0.008± 0.016 

0.00±0.00±0.03 
0.011±0.031 ±0.019 



Table 2. Results of /2, XE and XB ctt = 1 GeV^ for proton, neutron and p — n. The 
uncertainties are first statistical, then systematic. 



Target 


/2 


Xe 


Xb 


P 
n 

p — n 


-0.160 ±0.028 ±0.109 
0.034 ±0.005 ± 0.043 
-0.136 ±0.102 ±0.039 


-0.082 ±0.016 ±0.071 
0.031 ±0.005± 0.028 
-0.100 ±0.068 ±0.028 


0.056 ±0.008 ±0.036 
-0.003 ± 0.004 ± 0.014 
0.036 ±0.034 ±0.017 



The Fi moments were fit, varying the minimum Q^— value to study the conver- 
gence of the OPE series. The extracted quantities have large uncertainties (dom- 
inated by the systematic uncertainty) but are stable with respect to a minimal 
Q^— value when it is below 1 GeV^. The results do not vary significantly when the 
/i8 term is added, which justifies a posteriori the use of the truncated OPE series in 
the chosen Q^— range. In the proton case, the elastic contribution makes a signifi- 
cant contribution to the fie term at low but this does not invalidate a priori the 
validity of the series since the elastic contributes mainly to /ig and /is, but remains 
small compared to fi4. We notice the alternation of signs between the coefficients. 
This leads to a partial suppression of the h igher -twist effects and may be a reason 
for quark-hadron duality in the spin secto r^ ^^ l We also note that the sign alter- 
nation is opposite for the proton and neutron. Following Eq. (|73p . the electric and 
magnetic color polarizabilities were determined. Overall, the values given in Table 2 
are small, and we observe a sign change in the electric color polarizability between 
the proton and the neutron. We also expect a sign change in the color magnetic 
polarizability. However, with the large uncertainty and the small negative value of 
the neutron xb, it is difficult to confirm this expectation. 
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5.8. Quark-hadron duality in spin structure functions 

Detailed studies of duality in the spin structure functions gi and gf have been 
pubhshed by the CLAS egl collaboration. One observes a clear trend of strong, 
resonant deviations from the scaling curve at lower towards a pretty good agree- 
ment at intermediate Q^. The integral of gi over the whole resonance region begins 
to agree with the NLO results above « 1.7 GeV^ (see Fig. [71). The results on the 
proton and deuteron from eglb^^thus indicate a much slower approach to "global" 
duality for the polarized str uctur e function gi than has been observed for the un- 
polarized structure functions ^^^^ Local duality seems violated in the A— resonance 
region even for Q^— values as high as 5 GeV^. 

The data taken by the RSS collaboration in Hall C corroborate these obser- 
vations and add more precise data points for w 1.3 GeV^. 




13 3 4 5 r , , I , , , , I , , , , 

1 2 3 

(GeVVc^) Q-(GeV/cr 



Fig. 7. Data on gi{x,Q'^) averaged over the resonance region for the proton and the deuteron from 
Hall B egl (left panel). The hatched curves represent the range of extrapolated DIS results from 
modern NLO fits (GRSV and AAC), evolved to the of the data and corrected for target mass 
effects. The open circles in the left panel include the elastic contribution, while the filled circles 
are only integrated over W > 1.08 GeV. The right panel shows rf^'= and T'l of the resonance 
region from Hall A EOl-012, together with lower results from E94-010, compared with the 
world DIS fits. 



The spin structure functions gi and Ai were measured in the resonance region 
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(W^ < 2 GeV) for ^He in the -region below 1 GeV^ by Hall A experiments 
E94-010llJlISland from 1 GeV^ to 4 GeV^ by E01-012l^. Due to the prominent 
contributions from the A resonance, local duality does not appear to work at low 
(below 2 GeV^). At high (above 2 GeV^), the A- resonance contribution starts 
to diminish. The resonance data were integrated to study global duality. Figure [7] 
(right panel) shows the results for both "^He and the neutron in comparison with the 
DIS fits evolved to the same Q^. The resonance data agree with the DIS fits at least 
for higher than 1.8 GeV^, indicating that global duality holds for the neutron 
and ^He spin structure function, gi, in the high Q^— region (above 1.8 GeV^). 

The study of quark-hadron duality helps us to study higher-twist effects and 
decide where to extend the kinematic region apply partonic interpretations. From 
a practical point of view, the good understanding of the higher twist effects and 
quark-hadron duality allows us to considerably e xtend the experimental database 
used to extract the polarized parton distributions E^. 



5.9. The effective strong coupling at large distance 

In QCD, the magnitude of the strong force is given by the running coupling constant 
Us- At large Q^, in the pQCD domain, is well defined and can be experimentally 
extracted, e.g. using the Bjorken sum rule, (see eq. HI]). The pQCD definition leads 
to an infinite coupling at large distances, where approaches Ag^^^, This is not 
a conceptual problem because we are out of the validity domain of pQCD. Since 
the data show no sign of discontinuity when crossing the intermediate Q^— domain, 
see e.g. Fig. [3l it is natural to look for a definition of an effective coupling a'^J^ , 
which works at any and matches ag at large but stays finite at smal l Q^ . 
The Bjorken sum rule can be used advantageously to define al^^ at low 1113| 
The data on the Bjorken sum are used to experimentally extract a'^^^ following a 
prescription by Grunberg HH, see Fig. [51 The Bjorken and GDH sum rules also 
allow us to determine a'^J^ at, respectively, large and ~ 0. The extracted 
a'gf^ provides for the first time an experimental "effective coupling" at all . An 
interesting feature is that a^^-^ becomes scale invariant at small Q^, which was 
predicted by a number of calculations. 



6. Conclusion 

A large body of nucleon spin-dependent cross-section and asymmetry data have 
been collected at low to moderate in the resonance region. These data have 
been used to evaluate the Q^— dependence of moments of the nucleon spin structure 
functions gi and 52, including the GDH integral, the Bjorken sum, the BC sum and 
the spin polarizabilities, and to extract higher-twist contributions. The latter have 
provided access to the color polarizabilities in the nucleon. 

At close to zero, available next-to-leading order xPT calculations were tested 
against the data and found to be in reasonable agreement for — 0.1 GeV^ for 
the GDH integral /(Q^) and Ti{Q'^). Above = 0.1 GeV^ a significant difference 
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Fig. 8. Value of rxj^ j-K extracted from the data on the Bjorken sum, on r decay and on the 
Gross-Llewellyn-Smith sum. The values of ct%^ ^ computed within pQCD and using the Bjorken 
sum are given by the gray band. The values of ofj^ /n extracted using the GDH sum rule axe 
given by the red dashed line. 



between the calculation and the data is observed, pointing to the limit of appli- 
cability of xPT as becomes larger. The comparison of results for the forward 
spin polarizabilities 7o((3^) and 5lt with the xPT calculations show significant 
disagreements. 70 (Q^) has significant contributions from the A resonance and its 
inclusion in xPT calculations have not been fully worked out. The disagreements 
may indicate a better treatment of the A resonance is needed in the xPT calcu- 
lations. On the other hand, the xPT calculation of 5lt was expected to offer a 
faster convergence because of the absence of the A contribution. However, none 
of the available xPT calculations can reproduce 5lt for the neutron at the low- 
est point of 0.1 GeV^. This discrepancy presents a significant challenge to our 
theoretical understanding at its present level of approximations. 

Overall, the trend of the data is well described by phenomenological models. The 
dramatic Q^— evolution of Igdh from high to low was observed as predicted 
by these models for both the proton and the neutron. This behavior is mainly 
determined by the relative strength and sign of the A resonance compared to that 
of higher-energy resonances and deep-inelastic processes. This also shows that the 
current level of phenomenological understanding of the resonance spin structure 
using these moments as observables is reasonable. 

The neutron BC sum rule is observed to be verified within experimental uncer- 
tainties in the low-Q^ region due to a cancellation between the inelastic and the 
elastic contributions. The BC sum rule is expected to be valid at all . This test 
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validates the assumptions going into the BC sum rule, which provides confidence 
in sum rules with similar assumptions. 

In the Q^-region above 0.5 GeV^, the first moment of g-^ for the proton, neutron 
and the proton-neutron difference were re-evaluated using the world data and the 
same extrapolation method of the unmeasured regions for consistency. Then, in the 
framework of the OPE, the higher-twist /2 and d2 matrix elements were extracted. 
The low-Q^ data have allowed us to gauge the convergence of the expansion used 
in this analysis. The extracted higher-twist (twist-4 and above) effects are not 
significant for above 1 GeV^. This fact may be related to the observation that 
quark-hadron duality works reasonably well at above 1 GeV^ . 

Finally, the proton and neutron electric and magnetic color polarizabilities were 
determined by combining the twist-4 matrix element /2 and the twist-3 matrix 
element ^2 from the world data. Our findings show a small and slightly positive 
value of xe and a value of xb close to zero for the neutron, while the proton has 
slightly larger values for both xe and xb but with opposite signs. 

Overall, the new JLab data have provided valuable information on the transition 
between the non-perturbative to the perturbative regime of QCD. They form a 
precise data set for twist expansion a nalysis an d a check of xPT calculations. 

Recently completed measurements ^^^^^^^^ and future 12 GeV experiments 
on the (/I, (72 structure functions and d2 moments for the proton and neutron at 

« 1—5 GeV^ will reduce the uncertainty in the extracted higher-twist coefficients 
and provide a benchmark test of Lattice QCD. 

Acknowledgments 

Thanks to Alexandre Deur, Karl Slifer and Patricia Solvignon for providing fig- 
ures. Thanks to Alexandre Deur, Kees de Jager and Vince Sulkosky for careful 
proof-reading. This work was supported by the U.S. Department of Energy (DOE). 
The Southeastern Universities Research Association operates the Thomas Jefferson 
National Accelerator Facility for the DOE under contract DE-AC05-84ER40150, 
Modification No. 175. 

References 

1. V. Bernard, Prog. Part. Nucl. Phys. 60, 82 (2008). 

2. M. Gockeler et at, PoS LAT2006, 120 (2006). 

3. I. C. Cloet and C. D. Roberts, PoS LC2008, 047 (2008). 

4. See, e.g., J. Polchinski and M. J. Strassler, Phys. Rev. Lett. 88, 031601 (2002); ibid. 
S. J. Brodsky and G. F. de Teramond, Phys. Rev. Lett. 96, 201601 (2006); 94, 201601 
(2005). 

5. E. V. Shuryak and A. I. Vainshtein, Nucl. Phys. B 201, 141 (1982). 

6. E. Stein, P. Gornicki, L. Mankiewicz and A. Schafer, Phys. Lett. B 353, 107 (1995). 

7. X. Ji, Nucl. Phys. B 402, 217 (1993). 

8. O. Stern and W. Gerlach, Z. Physik 7, 249 (1921); ibid., 8, 110 (1921); ibid., 9, 349 
and 353 (1922). 



January 22, 2010 19:17 WSPC/INSTRUCTION FILE moments-DNPOS 



28 J. -P. Chen 

9. R. Frisch and O. Stern, Z. Physik 85, 4 (1933); I. Estermann and O. Stem, ihid., 7 
(1933). 

10. S. B. Gerasimov, Sov. J. Nucl. Phys. 2, 598 (1965) 

11. S. D. Drell and A. C. Hearn, Phys. Rev. Lett. 16, 908 (1966). 

12. R. Hofstadter, Rev. Mod. Phys. 28, 214 (1956). 

13. See, e.g., S. E. Kuhn, J. -P. Chen and E. Leader, Prog. Part. Nud. Phys., 63, 1 (2009). 

14. J. D. Bjorken, Phys. Rev. 148, 1467 (1966). 

15. JLab E97-110, J. P. Chen, A. Deur, F. Garibaldi, spokespersons; V. Sulkosky, Proc. 
Spin Structure at Long Distance, Edited by J. P. Chen, W. Melnitchouk and K. SHfer, 
AIP 1151, 93 (2009). 

16. M. Amarian et al, Phys. Rev. Lett. 89, 242301 (2002). 

17. M. Amarian et al, Phys. Rev. Lett. 92, 022301 (2004). 

18. M. Amarian et al, Phys. Rev. Lett. 93, 152301 (2004). 

19. K. SUfer et al, Phys. Rev. Lett. 101, 022303 (2008). 

20. P. Solvignon et al, Phys. Rev. Lett. 101, 182502 (2008). 

21. P. Solvignon, Proc. Spin Structure at Long Distance, Edited by J. P. Chen, W. Mel- 
nitchouk and K. Slifer, AIP 1151, 101 (2009). 

22. R. Fatemi et al, Phys. Rev. Lett. 91, 222002 (2003). 

23. J. Yun et al, Phys. Rev. C 67, 055204 (2003). 

24. K. V. Dharmawardane et al, Phys. Lett. B 641, 11 (2006). 

25. Y. Prok et al, Phys. Lett. B 672, 12 (2009). 

26. P. Posted et al, Phys. Rev. C 75, 035203 (2007). 

27. F. R. Wesselmann et al, Phys. Rev. Lett. 98, 132003 (2007). 

28. K. Slifer et al, arXivnucl-ex/08 12003 L. 

29. Z.-E. Meziani et al, Phys. Lett. B 613, 148 (2005). 

30. M. Osipenko et al, Phys. Lett. B 609, 259 (2005). 

31. A. Deur, arXiv nucl-ex/0508022 

32. A. Deur et al, Phys. Rev. Lett. 93, 212001 (2004); A. Deur et al., Phys. Rev. D 78, 
032001 (2008). 

33. J. Ahrens et al. Phys. Rev. Lett. 87, 022003 (2001). 

34. H. Dutz et al., Phys. Rev. Lett. 91, 192001 (2003). 

35. H. Dutz et al., Phys. Rev. Lett. 93, 032003 (2004). 

36. M. Anselmino, A. Efremov and E. Leader, Phys. Rep. 261, 1 (1995); ibid., 281 399 
(erratum) (1997). 

37. J. Soffer, O. V. Teryaev, arXiv | hep-ph/9906455| 

38. X. Artru, M. Elchikh, J. M. Richard7J. Soffer and O. V. Teryaev, ' hep-ph/0802016^ . 

39. D. Drechsel and T. Walcher, Rev. Mod. Phys. 80, 731 (2008). 

40. L. N. Hand, Phys. Rev. 129, 1834 (1963). 

41. P. Hoodbhoy, R. L. Jaffe and A. Manohar, Nucl. Phys. B 312, 571 (1989). 

42. L. L. Frankfurt and M. I. Strikman, Nucl. Phys. A 405, 557 (1983). 

43. D. Drechsel, S.S. Kamalov and L. Tiator, Phys. Rev. D 63, 114010 (2001). 

44. X. Ji and J. Osborne, J. of Phys. G 27, 127 (2001). 

45. M. Anselmino, B. L. loffe, and E. Leader, Sov. J. Nucl. Phys. 49, 136 (1989). 

46. D. Drechsel, B. Pasquini and M. Vanderhaeghen, Phys. Rep. 378, 99 (2003). 

47. D. Drechsel and L. Tiator, Ann. Rev. Nucl. Part. Sci. 54, 69 (2004). 

48. J. D. Bjorken and S. D. DreU, "Relativistic Quantum Fields", McGraw HiU, New 
York (1965). 

49. F. E. Low, Phys. Rev. 96, 1428 (1954). 

50. K. Wilson, Phys. Rev. 179, 1499 (1969). 

51. H. Burkhardt and W. N. Cottingham, Ann. Phys. (N.Y.) 56, 453 (1970). 



January 22, 2010 19:17 WSPC/INSTRUCTION FILE moments-DNPOS 



Moments of Spin Structure Functions: Sum Rules and Polarizabilities 29 

52. J. Ellis and R. L. Jaffe, Phys. Rev. D 9, 1444 (1974); ibtd., D 10, 1669 (1974). 

53. X. Ji and P. Unrau, Phys. Lett. B 333, 228 (1994). 

54. X. Ji and W. Melnitchouk, Phys. Rev. D 56,1 (1997). 

55. S. Wandzura and F. Wilczek, Phys. Lett. B 72, 195 (1977). 

56. M. Burkardt, Proc. Spin Structure at Long Distance, Edited by J. P. Chen, W. Mel- 
nitchouk and K. Slifer, AIP 1151, 26 (2009). 

57. J. Ashman et at, Phys. Lett. B 206, 364 (1988). 

58. SMC collaboration: D. Adams et ai, Phys. Lett. B 329, 399 (1994). 

59. SMC collaboration: D. Adams et ai, Phys. Lett. B 357, 248 (1995). 

60. SMC collaboration: D. Adams et al, Phys. Lett.B 396, 338 (1997). 

61. SMC collaboration: D. Adeva et ai, Phys. Lett. B 302, 533 (1993). 

62. SMC collaboration: D. Adeva et al, Phys. Lett.B 412, 414 (1997). 

63. K. Ackerstaff et al, Phys. Lett. B 404, 383 (1997). 

64. K. Ackerstaff et al, Phys. Lett. B 444, 531 (1998). 

65. A. Airapetian, et al, Phys. Lett. B 442, 484 (1998). 

66. M. J. Alguard et al, Phys. Rev. Lett. 37, 1261 (1976). 

67. M. J. Alguard et ai, Phys. Rev. Lett. 41, 70 (1976). 

68. P. L. Anthony et al, Phys. Rev. D 54, 6620 (1996). 

69. K. Abe et al, Phys. Rev. D 58,112003 (1998). 

70. K. Abe et al, Phys. Rev. Lett. 79, 26 (1997). 

71. K. Abe et al, Phys. Lett. B 404, 377 (1997). 

72. P. L. Anthony, et al, Phys. Lett. B 493, 19 (2000). 

73. P. L. Anthony, et al, Phys. Lett. B 458, 529 (1999). 

74. P. L. Anthony et al, Phys. Lett. B 553, 18 (2003). 

75. C. Bernet, Proceedings of DIS2005 (2005), Madison, Wisconsin, AIP Conf. Proc. 
(2005). 

76. A. Deshpande, Proceedings of DIS2005 (2005), Madison, Wisconsin, AIP Conf. Proc. 
(2005). 

77. X. Ji, Phys. Rev. Lett. 78, 610 (1997). 

78. M. Burkardt, Phys. Rev. D 72, 094020 (2005). 

79. HaU A collaboration: J. Alcorn et al, Nucl. Inst. Meth. A522, 294 (2004). 

80. HaU A Status Report - 2007. 

81. CLAS collaboration: B. A. Mecking et al, Nucl. Inst. Meth. A503, 513 (2003). 

82. L.W. Mo and Y.S. Tsai, Rev. Mod. Phys. 41, 205 (1969). 

83. I.V Akushevich and N.M. Shumeiko, J. Phys. G 20, 513 (1994). 

84. C. Ciofi degh Atti and S. Scopetta, Phys. Lett. B 404, 223 (1997). 

85. A. Abragam and M. Goldman, Rep. Prog. Phys. 41 396 (1978). 

86. D.G. Crabb and W. Meyer, Annu. Rev. Nucl. Part. Sci. 47, 67 (1997). 

87. C. D. Keith et al, Nucl. Inst. Meth. A 501, 327 (2003). 

88. L.W. Whitlow et al, Phys. Le tt. B 250, 193 (1990); 

89. Y. Liang et al., nucl-ex/0410027 

90. K. Abe et al, Phys. Rev. D 58 112003 (1998). 

91. M. Arneodo et al, Phys. Lett. B 364, 107 (1995). 

92. V. Burkert and Z. Li, Phys. Rev D 47, 46 (1993). 

93. J. Soffer, Phys. Rev. Lett. 74, 1292 (1995). 

94. N. Bianchi and E. Thomas, Nucl. Phys. B 82 (Proc. SuppL), 256 (2000). 

95. X. Ji, C. Kao, and J. Osborne, Phys. Lett. B 472, 1 (2000). 

96. V. Bernard, T. Hemmert and Ulf-G. Meissner, Phys. Lett. B 545, 105 (2002) 

97. V. Bernard, T. Hemmert and Ulf-G. Meissner, Phys. Rev. D 67, 076008 (2003). 

98. H. Krebs, et al, Proc. Spin Structure at Long Distance, Edited by J. P. Chen, W. 



January 22, 2010 19:17 WSPC/INSTRUCTION FILE moments-DNPOS 



30 J. -P. Chen 

Melnitchouk and K. Slifer, AIP 1151, 42 (2009). 

99. C. W. Kao, Proc. 18th Int. Spin Phys. Symp. Edited by D. Crabb, et al, AIP 1149, 
289 (2008). 

100. J. Soffer and O. V. Teryaev, Phys. Rev. D 70, 116004 (2004). 

101. V. D. Burkert and B. L. loffe, Phys. Lett. B 296, 223 (1992). 

10 2. The Science Driving the 12 GeV Upgrade of CEBAF, 
http: / /www.jlab.org/div_dept /physics_division/GeV.html. 

103. P. Mergell, Ulf-G. Meissner and D. Drechsel, Nucl. Phys. A 596, 367 (1996). 

104. S. A. Kulagin and W. Melnitchouk (to be published), private communication. 

105. C. W. Kao, T. Spitzenberg and M. Vanderhaeghen, Phys. Rev. D 67, 016001 (2003). 

106. M. Gockeler et al, Phys. Rev. D 63, 074506, (2001). 

107. V. D. Burkert, Phys. Rev. D 63, 097904 (2001). 

108. J. Blumlein and H. Boettcher, Nucl. Phys. 636, 225 (2002). 

109. K. Kramer et al, Phys. Rev. Lett. 95, 142002 (2005). 

110. X. Zheng et al, Phys. Rev. C 70, 065207 (2004). 

111. W. Melnitchouk, R. Ent and C. Keppel, Phys. Rept. 406, 127 (2005). 

112. E. Leader, A. V. Sidorov, D. B. Stamenov, Phys. Rev. D 63, 034023 (2006). 

113. A. Deur, V. Burkert, J. P. Chen and W. Korsch, Phys. Lett. B 665, 349 (2008);B 
650, 244 (2007). 

114. G. Grunberg, Phys. Lett. B 95, 70 (1980); Phys. Rev. D 29, 2315 (1984); Phys. 
Rev. D 40, 680 (1989). 

115. B. Sawatzky et al, Proc. Spin Structure at Long Distance, Edited by J. P. Chen, 
W. Melnitchouk and K. Slifer, AIP 1151, 145 (2009). 

116. O. A. Randon, Proc. Spin Structure at Long Distance, Edited by J. P. Chen, W. 
Melnitchouk and K. Slifer, AIP 1151, 82 (2009). 



